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We discuss the effect of the renormalization procedure in the computation of the 
unification point for running coupling constants. We explore the effects of threshold- 
crossing on the /3-functions. We compute the running of the coupling constants of 
the Standard Model, between mz and Mp, using a mass dependent subtraction 



procedure, and then compare the results with MS, and with the 9— function approx- 
imation. We also do this for the Minimal Supersymmetric extension of the Standard 
Model. In the latter, the bounds on susy masses that one obtains by requiring per- 
turbative unification are dependent, to some extent, on the procedure. 
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The last years have seen a revival of the study of scenarios for a possible perturbative 
unification of the strong, electromagnetic and weak interactions into a simple group G, the 
simplest candidate being SU(5) (5,0- In this way, the three couplings of the Standard 
Model (SM), g3,g2,g', for SU(3) C x SU{2)l x U(1)y, would evolve up to a common value g 
at the unification scale Mx (order of 10 16 GeV) from their disparate value at rriz- Thus the 
study of the evolution of these couplings from low to high energies, can be used to obtain 
information about the possible validity of the perturbative unification picture. 

Since we pretend to compare experimental values with theoretical predictions for the 
coupling constants gi, how to actually study the scale evolution of the couplings becomes 
an issue. The tool to carry this out is the renormalization group (RG), with /^-functions 
defined as /3j = from which we can get <7i(/i), /i being the energy scale at which 

we probe the system, and gi the renormalized coupling. The physical coupling constants 
are measured at the laboratory at a given energy scale /io, typically of order of mz- We 
identify these values with g^o), and with the help of the RGE we infer the evolution of 
the couplings as the scale changes. However, in order to compute the /^-functions, one 
needs to specify a renormalization scheme, and beta functions differ for one scheme to 
another. Broadly speaking, one has generally two choices of scheme at his/her disposal: the 
modified minimal subtraction scheme (MS) ||, and mass dependent subtraction procedures 
(MDSP) Q. In MS schemes, the /^-functions depend on the particle content at the energy 
scale on which one computes them and not explicitly on the masses nor energy scale. On 
the other hand, MDSP schemes take into account a dependence on particle masses: each 
graph including particles with mass rrii running in the loop, contains a function of the ratio 
fi 2 /m 2 associated with it, where \i is the energy scale. The decoupling theorem || ensures 
that the contributions from graphs with heavy field loops are suppressed at low momenta 
(much smaller than the masses involved). Naturally, the functions f(fi 2 /m, 2 ) associated 
with each graph have the property that f(fi 2 /m 2 ) — > 0, when ^ 2 /m 2 — > 0, and therefore the 
decoupling of the heavy degrees of freedom is explicitly enforced. Because of their origin, 
these functions include the threshold effects arising as \i becomes larger than the putative 
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particle mass. 

In addition to the above schemes, one can approximate the threshold effects by Heaviside- 
9 functions and impose this on the MS /3-functions, which are then integrated to give the 
effective couplings. Hence a massive particle only contributes at scales larger than its mass 
||. However, this is only an approximation to the contribution given by the functions 
/(/i 2 /m 2 ) in the MDSP, and the results for the effective gauge couplings can be appreciably 
different (the more thresholds are crossed, the more different are the effective couplings). 

To illustrate the relevance of the scheme-choice we calculate a e (mw), integrating the 
/3-function from m e with a~ l {m e ) = 137.036, to mw = 81 GeV and assuming m t > mz, 
using the three methods, viz., (a) MS without any reference to the thresholds (which of 
course gives extremely inaccurate results as it is shown below), (b) ^-functions, and (c) 
MDSP. It is straightforward to see that: 

"7 1 \ms= 120-06 , "e" 1 \e= 128.54 , a" 1 \mdsp= 129.33 . 

Because thresholds modify the derivatives of the couplings (/3— functions) and we are 
integrating over many decades in momenta, the effect due to the decoupling of massive 
degrees of freedom and the subtraction procedure could be important, and in any case one 
cannot disregard the effects of the thresholds. 

A further uncertainty is introduced when making considerations relative to unification, 
since the choice of the unifying group, G, impacts on the renormalization of the couplings. 
Let us assume that this group breaks into the SM at scale Mx- Since the particle content 
of G is different from the SM we will have, in general, light fields (particles of the SM) and 
heavy fields (with masses of the order of the unification scale Mx)- At low energies these 
heavy fields decouple from the theory. To take into account this decoupling in a proper way 
one can, e.g., integrate out the heavy fields from the action [[!]||. Carrying this out, the 
simple unification condition g^Mx) = g(M x ), is modified into #,(//) ~ 2 = g(^)~ 2 — 
where the functions Aj depend on the masses of the heavy fields, and the interpolating scale 
H satisfies rrij <C \x <C Mx- Again the values of the masses are important. 
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We are interested in studying how the presence of different masses affects the running 
of the couplings (including threshold effects), as well as its consequences on a possible 
unification scenario for the SM. Therefore we will study here the SM with a mass dependent 
renormalization procedure (MDSP). The relevant masses for the problem are mz,mw ,nk 
and rn/j, the Higgs mass. We do not consider the remaining fermion masses because they are 
much smaller than mz, so that for scales fi > rriz we can regard all the fermions as being 
massless except for the top quark. 

In the symmetric phase, when /i 3> mz, the couplings are g 2 , g\, g' 2 of SU C (3) x SUl{2) x 
£/y(l); when the symmetry breaks down, the relevant couplings are instead <?§, g\ and e 2 , 
with the coupling g' 2 is given as a function of g\ and e 2 by 

- = --- (1) 

9' 2 e 2 g 2 ■ [L) 

In the laboratory we can measure g 2 , g\ and e 2 but not g' 2 . Therefore we will calculate the 
evolution of these three couplings and infer the evolution of g' 2 from ([!]). 

We need to compute the /^-functions for g 2 , g\ and e 2 in an MDSP. In general, a 
renormalized coupling constant, g 2 , is related to the bare constant g^ by g 2 = Zg^, where 
Z is a product of renormalization constants. A choice of subtraction procedure as well as 
an arbitrary subtraction point /i, are implicit in the computation of the Z's. In an MDSP 



there are two contributions to Z: the pole part, typical of the MS procedure, and finite 
contributions which depend on both the masses of the virtual degrees of freedom and the 
subtraction point. Since the /3-function is defined as 

_ , dln g 2 _ 2 d\n Z 
dm fi dm /i 

this dependence on the masses and the subtraction point carries over to the /3-functions. 
Now, for a general gauge theory which includes fermions and scalars, extra complications 
arise: there will be different kinds of vertices involving the coupling constants (fermion- 
boson and scalar-boson vertex in the case of abelian gauge theories, and also three and four 
boson vertices for non-abelian gauge theories). For each of these vertices, Z is given as a 
product of different renormalization constants Zf 



z 



ex.leg^ v 




(3) 



Here Z3 is the wave function renormalization constant (WFRC) of the relevant gauge boson; 
Z v is the vertex renormalization constant; Z ex j eg are the WFRC of the external legs. The 
Ward-Takahashi identities guarantee the existence of relations among the infinite parts of 
different Z t , in such a way that the infinite part of Z does not depend on the vertex we 
choose. But, as it is well known, due to the presence of masses, this is not true for the 
finite part: in each vertex the masses involved are different. Yet another problem arises 
due to the gauge dependence of the Zf, although each Zi is gauge dependent, their product, 
Z, cannot depend on the choice of gauge since g% must be gauge independent. Again, this 
is true for the gauge dependence of the infinite part, but it is not true in general for the 
finite contributions. Therefore, in order to get an acceptable g% in an MDSP we have to 
define Z in a universal (independent of the kind of vertex), and gauge-independent way. 
This is easiest for abelian couplings, since in this case one can extend the abelian Ward 
identity, Z ex j eg = Z v , to the finite parts, in such a way that gf = Z%g\, Z 3 being both 
gauge independent and universal. This definition of g% is equivalent to defining an effective 
coupling via ||: 



where TIb is the transverse bare vacuum polarization. 

But this no longer works in the broken phase of non-abelian theories. Because of tree 
level mixing of the neutral gauge bosons of SU(2) and U(l), a calculation of radiative 
corrections to the proper self energies, Il/^n and Ilf^, at 1-loop order, gives a contribution 
Tl^zA) w hich is (i) not purely transverse and (ii) it further mixes the self-energies. In an 
self-explanatory notation: 



where U L (q 2 ) is in general proportional to a mass squared. The longitudinal term only 
appears when the symmetry is broken, since IT Miy is purely transverse in an unbroken gauge 



rFV) = GT? 2 - flWnV) + <rn V) , 



(5) 
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theory. Then, if the longitudinal part Hz a contributes to the mass matrix of the (Z°, A) 
system, the photon would not be massless at 1-loop order! To maintain a massless photon 
without further field redefinitions one needs to eliminate the contribution from H^a- 

These problems affect the couplings e 2 and gr|> since they reflect the non-abelian nature 
of the parent gauge theory. Kennedy and Lynn [|K]] have shown how to solve them, and, at 
the same time, define a g 2 . appropriate to phenomenology. They split the vertex contribution 
for g\ into a finite process-dependent part, and a "universal" part r , consisting of an infinite 
term (dictated by the Ward identities) and a finite term to be determined^. As pointed out 
by these authors, all the gauge dependence is included in the universal part. Before breaking 
the SU(2)l x U{1)y symmetry, the g\ coupling is redefined into 

<?2 2 = <7l(l + <7lr') • (6) 

When the symmetry is broken we can write the electric charge in terms of g~2 2 and g' 2 as 

e 2 g 2 2 g 12 

The advantage of working with the g 2 is now clear, since all non-abelian effects are now 
dumped into r , and one may define effective charges associated with the bare tilde param- 
eters in exactly the same way as it was done in the abelian case0, i.e. 

1 1 i / 9 



9i(<l 2 ) 9t 2 

For the couplings g 2 , e 2 we get the effective couplings: 

1 1 ,j 



n J B (g 2 ) . (8) 



g 2 (q 2 ) g 2 



(W B (q 2 ) - 2T (q 2 )) . (9) 



1 Prom now on, the term "vertex" will refer to the contributions due to both the vertex itself and 
external legs. 

2 Notice that the IL3 refers only to the transverse part 
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At the same time, if we calculate the 1-loop corrections to the mass matrix of the neutral 
bosons, the non diagonal term will be equal to H Z a + rri^rT , and this term is now finite. 
The finite part of r' may be chosen so that: 

U L ZA + m 2 w T' = . (10) 

Thus (a) the photon remains massless (at least to 1-loop) and (b) we have available 
scheme-independent effective charges 3 . 

For the SU (3) coupling the equivalent of (||) reads: 

^i + ( ngV)-2r<'V,, , (ii) 

where H^\q 2 ) is the bare proper self-energy of the gluon, and T^(q 2 ) is the universal part 
of the gluon vertex. Since gluons are massless, the process-independent part of the vertex 

(3) 

consists of an infinite and gauge dependent term, while U B includes the contributions due 
to massive quarks. 

The above equations for the effective charges can be written in terms of the couplings at 
scale mz in the form: 

-57^7 = 37^7 + ^W) - %(m 2 z ) - 2r'(g 2 ) + 2r'(m|)) I (12) 
the combinations H B (q 2 ) — H B (m 2 z ), T' (q 2 ) — V (m 2 z ) that appear in flT2|) are now finite. Using 



dimensional regularization we calculate IT B ,r'; in the Landau gauge we get the following 
expressions: 

(4n) 2 T' = 3 -Q-\n^f) + F r (a w ) , (13) 



3 It may be shown that (||) remains gauge independent when coupled to the RGE for the gauge 
parameter a; now, since a = is a fixed point of this RGE, we will choose to work in the Landau 
gauge where the value of a stays at zero, and therefore Goldstone bosons and ghosts remain 
massless. 
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n A = n^ + < + + n A + + n A , (14) 



(4vr) 2 n A 



WW 
2-nA 



(47r) 2 n^ + = -2aH,F 2 (a w ,0) 
1 /2 



11 = Sgll WA + C e ll WZ + + Sgll A+ + S ll Z+ + ll h+ + 11 2+ + lly , (15) 

13 /2 



(^rK A = -{^^-\nf-) +Fl (a w ,0)] , 
(47r) 2 n^ = - (J - ln ^) + F^, a z ) 
(47r) 2 n^ = -a w F 2 (a w , a fe ) , 
(47r) 2 n^ + = -a^F 2 (0,0) , 
(47r) 2 n^ + = -a z F 2 (a z ,0) , 
(4 7 r) 2 nr + = l^-ln^)+F 3 (a„0) , 



12 \e 

(47r) 2 n£: =-(--ln=# + - 
^ 12 V 

1 ^ /2 



(4-) 2 nf = ^ E (f - ^ ^ + F f (a fl , a /2 )) 



/ 

(4-) 2 r (3) = 7f--ln-# + ^) ■ (16) 



4 Ve 3, 

no) = n w + nf , 

13 /2 , . 07 



v ; 99 2 \e " 2 78/ 
(4-) 2 nf = ^E('-ln^ + ^/(«/,«/) 



where we have introduced the following abbreviations 



^? 2 -2/i 2 i 2 

a, ; = — - , s 9 = sm W , c 9 = 1 - s e 



2 2 

- = - — 7 + ln 47r . 

£ n - 4 
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The self-energies above come from the diagrams of Figs. ([!]) and (0). The sliding mass scale 
is denoted by /i, and n is the dimension of the space-time. Notice that the effective charges 
do not depend on fi, but only on the masses and the external momenta. The Fi functions 
contain the threshold effects for particle production. The functions i — 1, 3, /, behave 

like Ai ma, when aj goes to zero, and they tend to a constant value when aj goes to infinity. 
On the other hand, Fiiflj) — > in both limits. Therefore, in the high energy limit, dj — > 0, 



we recover the MS'-expressions: 

1 1 m| 

where, for 3 generations of fermions and 1 scalar doublet 

11 1 Q 

(4vr) 2 /3 e = - , (4tt) 2 /3 2 = --| , (4vr) 2 /3 3 = -7 . (18) 

Equation ( |l~2"D allows us to study the evolution of the couplings from rriz to higher energies, 
and therefore to extract consequences concerning perturbative unification, taking into ac- 
count all the caveats and ambiguities discussed above. The input values we will use, at the 



scale mz, are 11 



wa 2 d w = 0.2329 ±0.0013 , 
a- 1 = 127.9 ± 0.3 , 

a 3 = 0.111 ±0.003 , 
m w = 80.6 ± 0.4 GeV , m z = 91.161 ± 0.031 GeV . 

Due to the embedding of the SM in SU(5), one normalizes the U(1)y coupling as |fj|: 

-l 3 i 3 / _ x _ t \ 
a x = -a = -[a e L - a 2 ) . 

Since the top and the Higgs have not yet been detected, we will take their masses as free 
parameters, whose lower experimental bounds will be set as : m t > 90 GeV, rrth > 48 GeV 
|T2]] . On the other hand, these masses can not be much higher than about 200 GeV as 



required by perturbative bounds ||13|| . This is the value we will adopt for both m t and m^; 



however, it can be seen that values between 100 — 1000 GeV give rise to the same qualitative 
behavior. 

In Fig. (|J) we have plotted a' 1 ^) calculated with the three 1-loop procedures previously 
discussed: mass dependent (MD) procedure (eq. (|I2"D), integration with 9 functions, and 
MS. Only the high energy region (10 12 — 10 19 GeV), where unification 4 could occur, is 
depicted. We see at once that the approximation we use for «f 1 and ctjT 1 make some 
difference although in this case very little: for example, a^-MD is a little larger than the 
other couplings, due to the fact that when thresholds are taken into consideration the top 
quark decouples at low energies. This effect propagates to high energies by the RGE and in 
this case the decoupling is smoother than in the 9 approximation. 

In the running of a^ 1 a larger number of massive particles participate, and therefore 
their decoupling will introduce larger corrections than for the other two couplings. This is 
readily seen in the figure, which also shows that a^-9 * s a little larger than a^-MS (same 
reason as for a^ 1 ), and they both are larger than a^-MD. We now have the contributions 
of and Z°, which in the mass dependent scheme take "longer" to decouple. This is the 
dominant decoupling effect in , even though rriz, mw < m t, m h, which is not operative in 
ct2 X -9 because we begin to integrate precisely at the scale of mz- That is, the thresholds 
effects tend to make a^ 1 "less asymptotically free". This was already pointed out in Refs. 



14| , [fLq| which show that, within the context of the minimal SU(5) model, the /3-function 
for the SU{2) coupling with thresholds is positive in the region around mw 

We have also calculated a" 1 at 2-loop order without thresholds. The 2-loop RGE's in 



the MS are well known fll6| ; solving these equations by an iterative technique 
obtains: 



a -V)^rW>+(^>n^, (19) 



4 Of course this is not the case for the SM. This plot is simply intended to illustrate how the 
decoupling affects the evolution of the coupling constants. 
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«r 1 (/^) (1) = «r 1 (-.) + ^yin^ , 

where i, j — 1, 2, 3; the coefficients 6«, bij for 3 generations of fermions and 1 scalar doublet 
are: 



/ 41/10 ^ 



-19/6 
-7 



^ 199/50 27/10 44/5 ^ 



V 



(20) 



9/10 35/6 12 
11/10 9/2 -26 

In this case we calculate directly a^ 1 , rather than a~ x , because working in MS is equivalent 
to working with the symmetric theory SU(3) x SU(2) x U(l). 

In Fig. ([D we can see that 2-loop effects are qualitatively similar to MD effects: they 
raise a^ 1 and decrease , even though in the case of a^ 1 the effect from the decoupling of 
the gauge bosons dominates over 2-loop effects. Thus, if we take into account both of them 
(2-loop and thresholds) by calculating a" 1 with an MDSP at 2-loop order, a^ 1 will be less 
asymptotically free that a^ 1 -MS at 1-loop. 

Short of doing the exact calculation we have studied this case with an appropriate ap- 
proximation for the threshold functions. We approximate a~ l{ ~ l \ b^ in Eq. (|T9|) by: 



«r 1 (A*) (1) = *r\mz) + y^Ef (f W M l(- p2=m |) -f (k \a k ) | ( „ p2=At2) ) 

l 47r J k 

k 



(21) 
(22) 
(23) 



I k 



where we have summed over all the particles in the SM, and 



f {k \a k ) = ln(l + c k a k ) - \nc k a k 



f [k \ak) 



(24) 
(25) 



dhi(-p 2 ) l + c k a k 

These functions and f'^ have been chosen so that their behavior in the limit a k — ► 



is: 
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f k \a k ) -- In c k a k , - 1 , (26) 

and both functions vanish when a k — > oo. Thus in the high energy limit (or for massless 
particles) we recover the a^-MS expressions; the heavy masses decouple at low energy 
more smoothly than in the decoupling model by a step function. This approximation is 
very good at 1-loop order: comparing the numerical results with those obtained from a^ 1 - 
MD, the differences are only of order 0.1%. We also expect this to be true at 2-loop order. 
Actually it is not necessary to modify the 2-loop coefficient 6^, since there are no appreciable 



differences if we take fejj as given by the MS prescription or as given by Eq . (|23|) . This may 
suggest that 2-loop thresholds are not relevant (only 1-loop thresholds are). In fact, 2-loop 



threshold effects are quite small since they are corrections to the 2-loop MS contribution 
which is, in turn, a correction to the 1-loop contribution. 

If we compare the 1-loop and 2-loop results, both within this approximation, we see 
(Fig. that "2-loop+thresholds" has the effect of lowering a^ 1 and raising a^ 1 . Therefore, 
combined effects of mass and order of the perturbation theory, tend to bring a^ 1 and a^ 1 
closer at high energies. Unfortunately this effect is not strong enough to unify the three 
couplings of the SM, but it is interesting to point out its existence. 

It is clear this effect will become the more relevant the larger the number of massive 
particles we have. So far we have only worked with the standard matter content of the 
SM, in which all the particles have been detected, and their masses measured, except for 
the top and the Higgs. The particle contents can be modified, for example by extending 
the SM and working in the so called Minimal Supersymmetric Standard Model (MSSM). 
Here more degrees of freedom (the susy partners etc.) begin to contribute between mz 
and the putative unification scale. Since these susy particles have not been detected, there 
are at most lower bounds available to their masses. We can make use of some relations 



between them which reduce the number of the necessary arbitrary mass parameters [JTT . 
The common mass for both squarks and sleptons, m , and the gaugino mass, m^, may 
be bounded by demanding that susy masses be in the range between 45 GeV and 1 TeV. 



12 



We also take mz < m^, m + ~ itlh < 1 TeV, and m t , rrih — 200 GeV, where m M is the 
higgsino mass. 

Previous works in MSSM indicate that an Mx ~ O(10 16 ) and a common Mgusy (or 
explicit susy masses) in the range from mz to O(lTeV) are compatible with experimental 



bounds on proton decay [TT| , |TB| , |ID|] . In these analyses, and except for the case of 1-loop MS, 
the values for a 3 (m^) larger than 0.111 are favoured to have an Ms US y i n this range. The 



latest data on a^imz) from LEP (see table 1) [20[ indicate an ot^imz) larger than this by 
12%. Thus, until we have a more precise determination of ctzimz) we will have to take the 
allowed range for a^imz) to be (0.108, 0.125). Systematically one sees that the higher the 
value of a^iynz) the higher Mx, but the lower Ms US y This systematic also happens when 



we use 2-loop instead of 1-loop RGE's [|19j. The trend is maintained using Mass Dependent 
RGE's, but now the values of the susy masses needed to unify the couplings are higher than 
with the other methods. 

In Fig. ([]) we have represented a^-MD for two different values of the susy parameters 
and ctzivriz) = 0.111. We have also included a^-MS with Ms usy = rnz and the same value 
of as(mz) in order to compare the results yielded by different approximations. As pointed 
out before, the aJ l -MS unify, while the a^-MD do not. For example, we would need 
susy masses of order 10 2 TeV for a^rriz) = 0.111, and of order 1 TeV for a^imz) = 0.120. 
With a^(mz) in the range given before, and requiring Mx — O(10 16 ) to avoid conflicts with 
proton decay, susy masses of order m z are excluded, except m^j ~ which has to be 
less than 0(3 TeV). At 2-loop order, with complete light thresholds at 1-loop, the bounds 
increase: m 1 / 2 greater than approximately 3.5 TeV is excluded, but the rest of susy masses 
must be greater than O(10 TeV). 

So far, we have discussed perturbative unification within MSSM, but without any ref- 
erence to the unification group G, and the new heavy fields which are introduced by G in 
the theory. Moreover, we have required that the couplings "unify at the scale Mx" when 
they really "cut at the scale Mx" and run separately after it. To speak properly about per- 
turbative unification, we need that the couplings converge to only a single running coupling 
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constant: the G coupling. We would obtain this when we work within G, the group SU(5), 
taking into account the thresholds of all the masses, both light and heavy masses fli~4] . |21"f . 
As we have seen, light thresholds introduce appreciable differences in the running of the 
couplings, and the same effect takes place with the heavy masses when the scale approaches 
Mx [|22|| . Although we will have many more free parameters (the heavy masses), the trend is 
the same that in the simplest case when no reference was made to G: with mass dependent 
RGE's we need heavier susy masses in order to obtain unification with a^{mz) in the range 
allowed by experimental data. 

In conclusion, we have studied the evolution of the running couplings of the SM including 
complete threshold effects due to the light particles. We have also used effective charges fL0| , 



which are mass dependent, and process and gauge independent when we choose an adequate 
initial condition for the gauge parameter (Landau gauge). We have also seen that the effect 
of the decoupling of the light masses is not negligible at high energy. We apply the same 
method to study the MSSM, and the unification of the couplings within this model. Finally, 
We have found that light thresholds have the property of raising the values of the susy masses 
needed to keep a^{mz) within the range of the already available experimental values. 
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TABLES 



Experiment 


Central value 


Error 


ALEPH jets 


0.125 


±0.005 


DELPHI jets 


0.113 


±0.007 


L3 jets 


0.125 


±0.009 


OPAL jets 


0.122 


±0.006 


OPAL r 


0.123 


±0.007 


J/* 


0.108 


±0.005 


T 


0.109 


±0.005 


Deep Inelastic 


0.109 


±0.005 



TABLE I. Experimental values of a^(mz) 
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FIGURES 

FIG. 1. The 1-loop particle contributions to the H^) functions (j= A; W ± ). (a) Pure gauge 
boson contributions plus ghost (ui) contributions, (b) Graphs that mix scalars and gauge bosons, 
(c) Pure scalar contributions, (d) Fermion contributions. <p\ stands for the physical Higgs scalar; 
± , 4>2 for the Goldstone bosons associated with W ± and Z\ f for fermions (left and right); and 
fiL, Hl f° r the components of each doublet of left fermions. 

FIG. 2. The 1-loop particle contributions to the I^za function, i.e. to T' . (a) Gauge boson 
plus ghost, (b) Scalars and gauge bosons. 

FIG. 3. Evolution of the three couplings of the SM calculated with the three 1-loop procedures: 
mass dependent (solid lines), 6 function (dashed lines), and MS (dotted lines). 

FIG. 4. Evolution of the three couplings of the SM at 2-loop order with MS (solid lines), 
and 1-loop order with a mass dependent method (dashed lines). We also plot 1-loop MS for 
comparison (dotted lines). 

FIG. 5. Evolution of the three couplings of the SM at 1-loop order (solid lines), and 2-loop 
order (dashed lines), both with approximated threshold functions. 

FIG. 6. Evolution of the three couplings of the MSSM at 1-loop order calculated with (a) MS 
with M susy = mz (solid lines); (b) MD with m 1 / 2 = 45 GeV and mo = mz (dashed lines); (c) MD 
with mi/2 = mo = 1 TeV (dotted lines); we take = m + = mu = mo and m t = = 200 GeV. 
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